Derrick's theorem beyond a potential 
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Scalar field theories with derivative interactions are known to possess solitonic excitations, but 
such solitons are generally unsatisfactory because the effective theory fails precisely where nonlin- 
earities responsible for the solitons are important. A new class of theories possessing (internal) 
galilean invariance can in principle bypass this difficulty. Here, we show that these galileon theories 
do not possess stable solitonic solutions. As a by-product, we show that no stable solitons exist for 
a different class of derivatively coupled theories, describing for instance the infrared dynamics of 
superfluids, fluids, solids and some k-essence models. 



A famous theorem by Derrick states that in a theory of 
scalar fields with potential interactions only, there can- 
not be stable soliton solutions in more than one spatial 
dimension By 'soliton' it is usually meant a static 
non-trivial solution of the field equations with finite to- 
tal energy and in the present letter we stick to that def- 
inition. It is also understood that there is no external 
source sustaining the field profile - the solution should be 
sustained by non-linearities in the fields. Briefly speak- 
ing, Derrick's proof proceeds by contradiction: suppose 
a solitonic solution <po(x) exists, then consider deforma- 
tions 4>\(x) = 4>q(Xx) labeled by a dilation parameter A, 
and show that in more than one spatial dimension no 
(stable) stationary point of energy exists with respect to 
A for a scalar with purely potential interactions. 

What about interactions that are not purely of the 
potential form? After all, derivative interactions are 
generically expected in any low-energy effective theory, 
and beyond the two-derivative level one indeed intro- 
duces new powers of A into Derrick's argument, thus in 
principle allowing for more stationary points. The prob- 
lem, however, is that such interactions are always non- 
renormalizable. Therefore, either field gradients are so 
mild that they are negligible, or if they become impor- 
tant, the derivative expansion is expected to break down. 
In such a regime the effective theory is not enough to 
ascertain the existence of solitons, and one needs a UV 
completion (which comes with new or different degrees of 
freedom - thus in a sense redefining the problem we set 
out to solve) . For example, this is the case for skyrmions 
in the chiral Lagrangian [2j - to evade Derrick's theorem 
one needs to rely on terms beyond the two-derivative 
level, and there is no reason why further terms in the 
derivative expansion should be negligible. 

There is, however, a class of derivatively coupled the- 
ories where these conclusions do not necessarily apply. 
Consider the theory of a Goldstone boson it invariant 
under constant shifts of ir and its first derivatives 

7r(x) — > it + c + . (I) 

Such a symmetry has been dubbed 'galilean invariance', 



and 7r the 'galileon' 0. Galilean invariance forces the 
equations of motion to involve at least two derivatives 
acting on each field. Absence of ghosts in a non-linear 
regime demands that there be at most two derivatives on 
each field. Therefore, in order to have a galilean invariant 
dynamics that we can trust - at least classically - even 
when non-linearities are important, we need precisely two 
derivatives per tt in the field equations. This requirement 
corresponds to having Lagrangian terms of the form 



C n ~ <9tt <9tt (<9 2 7r) r 



(2) 



with suitable Lorentz contractions (and dimensionful co- 
efficients). Such combinations have been thoroughly clas- 
sified: in d + 1 spacetime dimensions there are d + 1, 
corresponding to n = 2, . . . , d + 2 0]. The n = 2 term 
is just the usual kinetic term (dir) 2 . That with n = 3 
is (<97r) 2 D7r, and so on. It is worth mentioning that the 
galileon describes the sub-horizon dynamics of the scalar 
sector of the DGP [3| model (n = 3), and of more generic 
theories that modify general relativity in the IR 0, S 0| ■ 
Of course classical consistency (e.g. absence of ghosts) 
of non-linear solutions is only a necessary condition 
for these to be inside the effective theory. One also 
needs quantum effects to be small. Here we are res- 
cued by galilean invariance, thanks to which the C n 's 
above do not get renormalized upon loop corrections, 
and terms with fewer derivatives per field are not gen- 
erated quantum- mechanically Q. Indeed the structure 
of divergences in the one-loop effective action (in 3 + 1 
dimensions) is schematically [6( 



r 1 - 1 -? ~ J2 [A 4 +A 2 a 2 +a 4 i og 9 2 /A 2 ] 



(3) 



where A is the mass-scale suppressing the galilean in- 
teractions ©■ For simplicity, we cut off the UV diver- 
gences at A. The sum runs over external legs. Now, 
a c i = ddir/A 3 is a measure of classical non-linearities. 
For instance the n-th order galilean interaction C n is 
roughly a™ ; -2 times the kinetic energy for n. On the 
other hand the quantity suppressing quantum effects is 
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really a q = <9 2 /A 2 , in the sense that, factoring out two 
powers of n, we have 



r i-ioop _ J2 [a q + a 2 + a 3 g loga q ]dTrdir 



ddir 



(4) 



Even for non-linearities of order one, ddir/A 3 ~ 1, this 
is suppressed w.r.t. to the tree-level interactions ^ by 
explicit powers of a q 14|. This means that for classical 



solutions with large non-linearities, quantum effects are 
small as long as gradients are mild d <C A 

These facts make it consistent (if not necessarily 'nat- 
ural') to postulate that the regime of large classical non- 
linearities is within the effective theory. The handful of 
galilean-invariant terms described above represent a con- 
sistent truncation of the theory in which we can self- 
consistently study non-linear classical solutions Q. It 
thus makes sense to ask whether there can be soliton 
solutions for the galileon. 

There is a final subtlety we need to address. As we have 
seen, in order to have small quantum effects we need mild 
gradients, d <C A. However, if the only scale appearing in 
the Lagrangian is A itself, and we introduce no external 
sources, it is clear that if a soliton exists it will have a 
size of order A -1 and gradients of order A. Nevertheless, 
we can consider the situation where different galilean in- 
teractions, say £3 and £4, are weighed by parametrically 
different scales, say A3 <C A4. If the soliton is sustained 
by a balance between these interactions, its size will be 
a see-saw combination of the two scales A3 1 and Aj 1 , 
and potentially parametrically bigger than either. On 
the other hand all our estimates above for r 1 "' 0013 ap- 
ply upon replacing A with the lower between the two 
scales, A3. In this example it is thus possible that there 
are classical soliton solutions that are within the effec- 
tive theory. Notice that the 'unnatural' tuning A3 A4 
is in fact technically natural, in that, as we mentioned, 
galilean interactions are not renormalized upon quantum 
corrections [15j . 

Let us now return to our main goal of discovering soli- 
tons, or lack thereof, in these theories. Suppose a soli- 
tonic solution ttq(x) exists. Small perturbations ip around 
it lead to energy fluctuations of this form: 
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d d xZ i: >{x)dnpdjip , 



(5) 



where Z 10 is a matrix built out of second derivatives of 
the background solution ir @, 0] , 



Z* j = c 2 S lJ + c 3 ((5^ V 2 7r - cWvro) + 



(6) 



where the c„'s are the coefficients of the terms ([2]) in the 
galileon action. This is highly non-trivial. In principle, 
given the structure of a generic galilean term @, one 
expects contributions to the quadratic action for small 



fluctuations of three different forms: 

57r 57r (a 2 7r )"- 4 d 2 ipd 2 ip 
dTT (d 2 TT Q ) n ' 3 d 2 ipd<p 
(<9 2 7ro)™~ 2 dtpdip 



(J) 
(8) 
(9) 



with suitable index contractions. Only the third struc- 
ture is of the form we like, eq. (0. However the second 
structure can always be rewritten as the third by integra- 
tion by parts, since the d 2 ipd(p piece is a total derivative: 

dftdv^datp = d^(d v )(pd a (p) ~ ±<9 Q (<9^<9„<^) . (10) 

As for the first structure, it cannot arise because the 
equations of motion are second order. Indeed for a 
non-trivial Lagrangian term A^ va ^ d^d^ip d a dp(p, there 
is an associated four-derivative contribution to the tp 
e.o.m. 2A^ a P d^d^dad/sip, which vanishes only if the 
tensor A does. Since the e.o.m. for the full galileon field 
are of second order by construction, a term of the form ([7]) 
cannot appear in the quadratic action for fluctuations. 

The quadratic structure in eq. ([5]) takes on signifi- 
cance in light of another fact, namely the existence of 
zero modes, i.e. small perturbations with SE = 0. The 
simplest example is the translational mode, a perturba- 
tion of the form: 



Vtt 



(11) 



where e is some infinitesimal vector. In other words, spa- 
tially displacing a soliton does not change its energy. 

Now, recall the quadratic structure for SE: for the 
soliton to be stable (or marginally so), it is crucial that 
Z %3 (x) be a positive semi-definite matrix everywhere in 
space. At large distances from the 'core' of our soliton, 
where ttq becomes very small, Z lJ is dominated by the 
zeroth order term, which is positive definite for ci > 
(eq. (J6J) ) . But the existence of a zero-energy perturbation 
then implies that in some other region Z 1 ^ must develop 
a negative eigenvalue: since the integrand in eq. §5§ is 
strictly positive at large x, somewhere else it should be- 
come negative for the integral to yield zero. 

The existence of negative eigenvalues for Z lJ in some 
region means that suitably chosen localized perturbations 
can be found that lower the energy of the solution - 
it is enough to pick very short wavelength wave-packets 
with momentum along the negative eigenvalue direction 
Such an instability, which is associated with neg- 
ative gradient energy for certain fluctuations, is much 
worse than an instability associated with a negative mass 
term. This is because the former also plagues very short- 
wavelength fluctuations, down to the UV cutoff of the 
theory. The decay rate is thus extremely fast, dominated 
by the shortest scales in the theory, and cannot be reli- 
ably computed within the effective theory. Conversely, a 
standard tachyon-like instability, like Jeans's or that as- 
sociated with the negative mode of a true-vacuum bub- 
ble, is dominated by modes with momenta of order of 
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the tachyonic mass scale, which can be parametrically 
smaller than the UV cutoff. As a consequence a tachyon- 
like instability can be slow, and its evolution can be con- 
sistently studied inside the effective theory, along with 
its interesting phenomenological consequences 16|. 

In conclusion, there is no consistent soliton solution 
in the galileon theory, in any number of spacetime di- 
mensions. Notice that our proof of instability crucially 
relies on the purely kinetic structure of the quadratic en- 
ergy for small fluctuations, eq. <j5j> - So, it does not apply 
to other theories that are known to possess soliton solu- 
tions, like e.g. a scalar theory with a potential in 1 + 1 
dimensions. There, the soliton has zero-energy transla- 
tional modes, yet the vanishing of their energy is accom- 
plished by having a localized negative mass term for small 
fluctuations. As a consequence, if one tries to construct 
negative-energy fluctuations by localizing them where the 
mass term is negative, one in fact enhances the positive- 
definite kinetic energy - which dominates over the mass 
term for large gradients - thus ending up with positive 
overall energy. Indeed one can show that for 'kinks' in 
1 + 1 dimensions the energy spectrum of small fluctua- 
tions in bounded by zero from below Q. 

As an interesting aside, it is worth pointing out that 
were the galileon soliton to exist, it would have zero total 
energy. To see this, consider the galileon action which 
takes the form S = J d d+1 x c„£ n (eq. ©). The 
total energy of a static configuration is 17 1 



E=- Jd d X Y,CnCn =J2E n , 



(12) 



and the soliton solution ttq(x) we are looking for should 
be a local minimum of this. Now we apply a minor gen- 
eralization of Derrick's argument. Consider the field con- 
figuration obtained by rescaling the x-dependence of ttq 
as well as its overall normalization 

= w7r (Ax). (13) 

A necessary condition for ttq to be a solution is that 

d x E(\,w)\ (11) =Q, d u E(X,u)\ (ltl) =Q, (14) 

where E(X,ui) is the energy of tt\ m . Now, the n-th order 
galilcan invariant term ([2]) involves n fields and 2n — 2 
derivatives. We thus have 



E n (X,u) = \ 2n - 2 - d uj n E^ 



where En is the n-th order term in the energy of 7Tq. 
Each E n and therefore their sum obeys 

Xd x E = -{d+2)E + 2ud u E . (16) 

We thus see that if a stationary point of E exists, it 
must have E — 0. This is not impossible a priori in our 
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X d+2 



(X^E^, (15) 



theory. Due to the peculiar higher-derivative structure 
of the action, the galileon can violate the null energy 
condition without obvious pathologies in the low-energy 
effective theory 0. This in principle allows for localized 
non-trivial field configurations with vanishing - even neg- 
ative - total energy. 

Incidentally, eq. (ITBl implies that there is another set 
of zero modes, namely deformations described by eq. fj 13[) 
but restricted to X 2 uj = 1. Such a deformation changes 
the overall energy by 1/X d+2 , but since the soliton (if 
it exists) has vanishing E, so does its deformation. In 
infinitesimal form (A = 1 + e), the deformation is (p e = 
— 2e 7To + e x- \7ttq. Thus our proof of the instability of the 
galileon soliton could have made use of this zero-mode 
instead of the translational zero-mode. 

Let us close with a discussion of another class of deriva- 
tively coupled theories for which essentially the same ar- 
guments apply, 



S= / d d+1 xP({d4>) 2 ) , 



(17) 



where P is an arbitrary function. Such an effective the- 
ory describes for instance the low-energy dynamics of a 
superfluid [Toj j . or, with more than one field, of ordinary 
fluids and solids [ll[ (our arguments also apply for mul- 
tifield generalizations of (fT7| ). It also describes more ex- 
otic systems, like the ghost condensate [HI], special cases 
of k-essence [13j , and simpler ones like a Goldstone boson 
non-linear ly realizing a U(l) symmetry, (j) — > <f> + c. As- 
suming there exists a soliton solution, one can see that 
Xuj = 1 is the flat direction in this case, on top of the 
translations ones of course. Provided there are spatial 
regions where small fluctuations carry a positive energy, 
as they should far away from the supposed soliton, there 
should also be regions where suitably chosen localized 
perturbations can destabilize the soliton. Indeed it is ob- 
vious that the perturbations' quadratic energy is still of 
the form {H} , of course now with a different Z l i . 

Notice that most of the systems mentioned above as 
examples for eq. (|17p spontaneously break Lorentz in- 
variance. So the physical question in this case is whether 
there are soliton solutions in the broken phase. Yet our 
proof never makes use of Lorentz invariance, so it applies 
unaltered in the broken phase. Perhaps more relevant 
is the worry for spontaneous breaking of time transla- 
tions; a superfluid, the ghost condensate, and k-essence 
all break time translations, thus making the definition of 
energy in the broken phase more subtle than usual. How- 
ever there is a linear combination of time translations and 
shift symmetry (0 — > <f> + c) that is unbroken. The cor- 
responding Noether charge is a perfectly good energy for 
excitations in the broken phase [12| , which our soliton so- 
lution should minimize, and which, apart from having a 
different 'tensor' structure, has the same scaling proper- 
ties as the original energy, thus lending itself to our proof. 
It should be emphasized however that the ghost conden- 
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sate is degenerate at the lowest derivative level, so that 
the excitations' gradient energy is in fact dominated by 
higher- derivative terms, of the form (V 2 7r) 2 , while the 
interactions are those given by eq. (fTTf and therefore 
have one derivative per field. This mismatch impairs the 
simple kinetic structure of the quadratic fluctuation La- 
grangian (0. Therefore our proof does not hold for the 
ghost condensate. 

Finally let us comment on the quantum properties of 
a theory like eq. ([T7|) . or of its multi-field generalizations 
describing fluids and solids. Consider the structure of 
divergences in the one-loop effective action. Assuming for 
simplicity that P is a generic function where all powers 
of (d(f>) 2 are weighed by the same scale A, we have 

r 1 -*** ~ ]T [a 4 + A 2 d 2 + d A log A 2 /d 2 ] (^P) m , 

(18) 

where we cut off the UV divergences at A. The quartic di- 
vergence renormalizes by order one terms already present 
in P((d(f>) 2 ). The other contributions - the quadratic di- 
vergence and the logarithmic one - involve more deriva- 
tives per field. That is, apart from a trivial renormal- 
ization of the tree-level Lagrangian, here too quantum 
effects are negligible even for solutions where (d<fi) 2 has 
large overall variations, as long as derivatives of (d(f) 2 
are everywhere small. The same is true at all loops as 
well. We even have experimental evidence that this is 
correct for certain systems - we can easily compress a 
fluid to a fraction of its original volume (A(dcf>) 2 ~ A 4 ) 
without exiting the effective theory, provided we do so 
slowly enough (d <C A). Like for the galileon, if A is 
really the only scale in the theory, our hypothetical soli- 
ton would have a size of order A -1 , and quantum effects 
would be large. We can bypass this problem by having 
very different scales in the theory and a see-saw mech- 
anism, as discussed above. But this requires a tuning, 
which here is not technically natural because, unlike for 
the galileon, the tree-level Lagrangian does get renormal- 
ized. Our argument shows that, even allowing for such a 
tuning, there cannot be stable solitons in these systems 

is]. 

Acknowledgements. The authors would like to thank 
S. Dubovsky, R. Rattazzi and especially E. Weinberg for 
very useful discussions and comments. LH thanks NYU 
and IAS for hospitality. This work is supported in part 
by the DOE (DE-FG02-92-ER40699) and NASA ATP 
(09-ATP09-0049). KH acknowledges support by funds 
provided by the University of Pennsylvania. 



[1] G. H. Derrick, J. Math. Phys. 5, 1252 (1964). 
[2] T. H. R. Skyrme, Nucl. Phys. 31, 556 (1962). 
[3] A. Nicolis, R. Rattazzi and E. Trincherini, 
larXiv:0811.2T97l [hep-th] . 



[4] G. R. Dvali, G. Gabadadz e and M. Porra ti, Phys. Lett. 
B 485, 208 (2000) [arXiv:hep-th/0005016] . 

[5] M. A. Lut y, M. Porrati and R. R attazzi, JHEP 0309, 
029 (2003) [arXivihep th/0303116] . 

[6] A. Nicolis and R. Ra ttazzi, JHEP 0406, 059 (2004) 
[arXiv:hep-th/0404159| . 

[7] A. Adams, N. Arkani-Hamed, S. Dubovsky, A. Nico- 
lis and R. Rattazz i, JHEP 0610, 014 (2006) 
[arXiv:hep-th/0602178] . 

[8] A! Nicolis, R Rattazzi and E. Trincherini, 
larXiv:0912.4258l [hep-th] . 

[9] S. R. Cole man, Subnucl. Ser. 13, 297 (1977). 
[10] D. T. Son, |arXiv:hep-ph/0204199l 

[11] S. Dubovsky, T. Gregoire, A. Nicolis and R. Rattazzi, 
JHEP 0603, 025 (2006) [arXiv:hep-th/0512260| . 

[12] N. Arkani-Hamed, H. C. Cheng, M. A. Luty and S. Muko - 
hyama, JHEP 0405, 074 (2004) [arXiv:hep-th/031 2099 . 

[13] C. Armendariz-Picon, V. F. Mukhanov and 
P. J. Steinhardt, Phys . Rev. Lett. 85, 4438 (2000) 
[arXiv:astro-ph/0004134l . 

[14] Strictly speaking, calculable (and measurable) quantum 
effects are only those associated with log-divergences - 
the last term in brackets. Power-divergences are instead 
completely regularization dependent, and indeed in di- 
mensional regularization they do not even show up. How- 
ever here we also want to stress that our galilean La- 
grangian does not even suffer from fine-tuning problems, 
which are associated with power divergences. 

[15] It is worth mentioning that a 'standard' UV-completion 
- i.e. a relativistic UV completion obeying the standard 
analyticity properties of S-matrix theory, like renormal- 
izable Lorentz-invariant QFTs and weakly coupled string 
theories - does not appear to exist for galilean theories 

an. 

[16] Having a flat direction - a trajectory in field space along 
which E remains constant - is not by itself necessarily a 
sign of instability. In the case of a flat direction in a multi- 
field potential, there is a continuum of different config- 
uration all with the same energy, but to go from one to 
the other, one should change the boundary conditions at 
infinity. Instead, local excitations have positive-definite 
energy coining from the field gradients. So, in that case, 
the flat direction does not correspond to localized per- 
turbations of the solution, and consequently each point 
along the flat direction corresponds to a stable solution. 
However in our case the situation is very different. The 
translation zero-modes, eq. (|11[1 . are local perturbations 
because no is localized. 

[17] The fact that for static configurations the Hamiltonian 
density is minus the Lagrangian density may seem non- 
trivial here, due to the higher-derivative nature of the 
action. In fact the Too one gets by deriving the action 
w.r.t. the metric is not —£.. On the other hand the Too 
derived a la Noether is always — C for static fields. The 
difference between the two is a total spatial derivative, 
which integrates to zero for localized configurations. 

[18] Localized, stationary vortices in a fluid do not conform 
to our definition of solitons: for them the velocity field is 
indeed time-indpendent - yet it involves a time deriva- 
tive of the canonical variables (f) 1 appearing in the La- 
grangian. That is, for a vortex the <^ 7 's have a non-trivial 
time-dependence [ill ]. 



